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A constituent parton picture of hadrons with logarithmic conhnement natu- 
^ ■ rally arises in weak coupling light-front QCD. Conhnement provides a mass 

gap that allows the constituent picture to emerge. The effective renormalized 
Hamiltonian is computed to 0{g^), and used to study charmonium and bot- 
tomonium. Radial and angular excitations can be used to fix the coupling 
CN \ a, the quark mass M, and the cutoff A. The resultant hyperfine structure is 

' very close to experiment. 
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The solution of Quantum Chromodynamics in the nonperturbative do- 
main remains one of the most important and interesting unsolved problems 
in physics. The basic assumption upon which our work is based is that it is 
possible to derive a constituent picture for hadrons from QCD If this is 

possible, nonperturbative bound state problems in QCD can be approximated 
as coupled, few-body Schrodinger equations. 

To arrive at a constituent approximation to QCD, we first want to separate 
vacuum fluctuations. This is achieved by formulating the theory on the light- 
front where the vacuum is trivial in presence of cutoffs. We start with the 
canonical light-front Hamiltonian, regulated by a cutoff on light-front energy, 
and use the similarity renormalization group to renormalize We expect the 
Hmailtonian to contain novel finite counterterms when there is spontaneous 
symmetry breaking [Q]. The vacuum is unchanged by spontaneous breaking. 

Renormalization of light-front Hamiltonians is more complicated than that 
of Lagrangians because many symmetries are not kinematically manifest, be- 
cause there is a separate power counting for longitudinal and transverse di- 
rections, and because locality is violated in the longitudinal direction. 

The basic idea of the similarity renormalization group is illustrated in 
Fig.l. Fig. 1(a) schematically shows a regulated bare light-front Hamiltonian 
matrix. The regulator used for illustration is a cutoff restricting light-front 
energies, and it makes the Hamiltonian matrix finite. There are different ways 
to impose the cutoff, but the specifics are not important for this discussion. 
The bare Hamiltonian contains couplings between all energy scales. This is 
the source of ultraviolet divergences. In order to renormalize the Hamiltonian, 
one needs to find counterterms that remove dependence on the cutoff. The 
similarity renormalization group is based on the following observation: if a 
Hamiltonian has a band diagonal form, as the Hamiltonian in fig. 1(b), then 
no ultraviolet divergence can appear at any finite order of perturbation theory 
because of the finite width of the Hamiltonian, as long as its matrix elements 
are finite. Therefore, if one can find a similarity transformation that trans- 
forms a bare Hamiltonian Hb as in fig. 1(a) to a band diagonal Hamiltonian 
H' as in fig. 1(b): 

H' = SHbS^ + counterterms, (1) 

it is possible to identify counterterms by requiring that the matrix elements of 
the band diagonal Hamiltonian are independent of the regulator. This ensures 
that physical observables will also be independent of the regulator. 

The similarity renormalization can be done in steps or continuously [^. 
In what follows, we use a discrete perturbative formulation around the free 
light-front Hamiltonian. 
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FIG. 1. (a) An example of a bare regulated Hamiltonian in light-front energy space. Energies 
run from zero up to an initial cutoff which is indicated by the subscripts of the first and the 
last diagonal elements, a in this figure schematically indicates nonzero matrix elements, but the 
matrix elements are not necessarily equal. The Hamiltonian couples states of all energy scales, (b) 
An example of a band diagonal Hamiltonian. The energies run up to the initial cutoff, but the 
Hamiltonian couples only states which are close in energy. 



With the divergences removed, the remaining task is to adjust finite parts 
of the counterterms. This is, in principle, achieved by restoring Lorentz in- 
variance and other exact symmetries in physical observables. However, if the 
similarity transformation can be done analytically, as in the calculation pre- 
sented here, it is straightforward to use coupling coherence which uniquely 
fixes all counterterms without explicit reference to underlying symmetries. 
The basic idea of coupling coherence is the following: in the Hamiltonian re- 
stricted by symmetries, albeit not manifest, the strengths of all operators are 
not independent but depend only on a finite number of independent canoni- 
cal parameters. Under a full renormalization group transformation (including 
change of scale and rescaling), the Hamiltonian reproduces itself in form ex- 
actly, apart from the change of the explicit cutoff and the running of those 
few independent couplings. All dependence on the cutoff is absorbed into the 
independent running couplings. Once one obtains a Hamiltonian that repro- 
duces itself as the cutoff is lowered and subtracts the divergences, any initial 
cutoff can be sent to infinity. 

The coupling coherent solution at second order in a generic interaction v 
(around the free light-front Hamiltonian ho) is 

Hab = {a\ho + v\h) 



J2k VakVkh 



■){\l^ak\-^)e{\^ak\-\^bk\) , e(|A 



-Aafc 



+ 



-)e(|A6fc|-[A,fc|) 



(2) 
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where Aij = Eqi — Eqj is the difference in light-front free energies, is the 
similarity scale, and sum over states k is limited by the initial cutoff and the 
explicit similarity cutoff in eqn. (2). A self-energy counterterm is also needed 



We would like to note that the similarity renormalization scheme which 
results in band diagonal effective Hamiltonians is a renormalization method 
that brings us closer to the desired constituent picture of hadrons. Indeed, if 
the constituents are massive, then at any finite order of perturbation theory 
with the band diagonal Hamiltonian, there is only a finite number of Fock 
states that couple to the lowest Fock component in the hadron. In our work 
low energy gluons acquire a mass gap due to a mass counterterm and due to 
the confining interaction, and high multiplicities of gluons are suppressed. 

After renormalization is completed, one is left with the effective Hamilto- 
nian band diagonal with the width of a hadronic scale. The effective Hamilto- 
nian contains complicated potentials, which result from eliminating the cou- 
pling between high and low energy states. It still contains emission and ab- 
sorption interactions, but these no longer mix states of high and low energies. 
However, if one tries to diagonalize the effective Hamiltonian directly, the 
wave function of a hadron must contain arbitrarily many parton components. 
Instead, we divide the effective Hamiltonian into a part Hq which is solved 
nonperturbatively, and the remaining part, V, is treated in bound state per- 
turbation theory. The division is arbitrary, but a choice of Hq missing an 
important part of physics of the system under consideration would lead to 
divergent bound state perturbation theory. Therefore, we want to choose Hq 
that approximates the physics relevant for hadronic bound states as closely 
as possible, and at the same time we want it to be manageable. We take a 
hint from the constituent quark model and include constituent masses and 
two-body potentials produced by the similarity transformation. 

A major simplification is achieved by not including emission and absorp- 
tion of low energy gluons in Hq. Once the particle- number-changing interac- 
tions are put in V, different Fock states decouple to leading order. The Hamil- 
tonian Hq provides an approximate qq valence quark description of mesons. 
The errors in approximation can be determined from bound-state perturba- 
tion theory in V which links the qq to multi-body Fock states. Mixing of 
different Fock components first enters at second order bound state perturba- 
tion theory. The valence approximation is best justified for heavy quarkonia. 

Based on the success of the constituent quark model, it is reasonable 
to choose a nonrelativistic limit of the effective Hamiltonian for Hq. This 
approximation, too, is best justified for heavy quarkonia. 

We have already used the effective Hamiltonian approach to study prop- 
erties of heavy-light mesons, in particular B mesons. In this letter we present 
numerical results obtained by applying the approach to charmonium and bot- 
tomonium for which the approximations are better justified. We fit IS, IP 
and 2S levels for both systems. We then predict hyperfine splitting in the 
charmonium ground state. The prediction is in good agreement with experi- 
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ment. 

We find the effective Hamiltonian to 0{g\). The effective Hamiltonian, 
which is generated by the similarity transformation and couphng coherence 
to order g"^, is band-diagonal in light-front energy with respect to a hadronic 
scale and it can be written as 0: 

^^cff = -f^free + ^1 + + t'2 eff , 

where -f^free is the light-front kinetic energy (we remind the reader that the 
light-front kinetic energy of a particle with transverse momentum and lon- 
gitudinal momentum p"*" is ^ ), vi is 0{g) emission and absorption with 
nonzero matrix elements only between states with energy difference smaller 
than the hadronic scale 

Let Pi, ki be the light-front three-momenta carried by a quark and an an- 
tiquark; Uj, Aj are their light-front helicities; u{p, a), v{k, A) are their spinors; 
index i = 1,2 refers to the initial and final states, respectively. Let q = pi—p2 
be the exchanged momentum and q~ = ^q^. V2 in eqn. (2) is an 0{g'^) 
instantaneous interaction with the following matrix element for free states 
containing a quark and an antiquark: 

- 9iu{p2: cr2h^u{pi,ai)v{k2, \2)Yv{ki, Xi){TaTb) 
X ^V^.Vu - \{pi + k{) - + A:2-)|) , 

where r/^w^ = defines the unit vector rj^. V2 off includes the 0{g'^) effective 
interactions generated by the similarity transformation. The effective inter- 
actions generated to this order contain one-body and two-body operators. In 
particular, the effective one-body quark operator is: 

"A^^ \ ,o^\2.„. -l^^' 



2— log -_ + 2— log 



2V+ ^2 a;„^A2 + M2 
M2 M2 1 



+3 log 



Xa Xa^A^ + M2 



where Xa = ^ is the longitudinal fraction of the momentum carried by the 
constituent under consideration, M is its mass, P"*" is the total longitudinal 
momentum of the state, is the longitudinal scale required in the cutoff 
by dimensional arguments, and e is an infrared cutoff which is to be taken 
to zero. The divergence in the effective one-body operator exactly cancels 
against a divergence in the effective two-body operator if the state is a color 
singlet §. 

The effective two-body operators have the following matrix elements be- 
tween states containing a quark pi and an antiquark: 



g\u{p2, (T2)l^u{pi,ai)v{k2, \2h''v{ki, Xi){TaTb) 



ei\Di\- 



+ 



-)e{\D2\-\D^\) 



D2 



where D^ij{q) 



±2 / 



g}[^ is the gluon propagator 



in hght-front gauge, Di, D2 are energy denominators: Di = Pi — P2 — 
and D2 = k2 — — ■ It has been shown that H^s contains a logarithmic 
confining interaction in addition to the Coulomb interaction [Q]. 

This is the output of the second order similarity transformation for qq 
matrix elements. 

For the purpose of bound state calculations, we split the effective Hamil- 
tonian (4) into Hq, which is solved nonperturbatively, and V = H^.^ — Hq. 
First, we make a nonrelativistic reduction of the effective Hamiltonian (4). 
In the nonrelativistic limit, the light-front scale ^ is naturally replaced by 
C = ^ ^ , where M is the mass of the heavy quark, and C carries the di- 
mension of mass Further, light-front momenta are naturally replaced by 
center-of-mass equal-time momenta in the nonrelativistic limit 

The spin-independent part of the two-body effective interactions includes 
a short-range Coulomb potential and a rotationally noninvariant long-range 
logarithmic potential. The confining potential arises due to an incomplete 
cancellation allowing the gluons becoming nonperturbative. The confining 
potential is a complicated function but in the nonrelativistic limit it can be 
double Fourier transformed (for longitudinal and transverse separation), and 
expanded in even Legendre polynomials |^. 

For Hq we choose the nonrelativistic reduction of the kinetic energy, the 
effective one-body operators, Coulomb potential and rotationally symmetric 
part of the confining potential with constituent masses. The Hamiltonian Hq 
is: 

1 



Hq = 4M 



2m 



\ Vo{Cr) 



TT 



where m is the reduced mass, r is an equal-time separation between the quark 
and the antiquark, and Vo(£r) is the angular average of the confining potential 
generated by the similarity transformation. It depends only on the separation 
of the quarks: 



Vo{Cr) = 21og7^ - 2Ci{n) + 4 



Sim (I 



n 



- cos IZ) sin TZ 



27 



where TZ = Cr and 7 is Euler constant. S contains the finite shift produced 
by the self-energies after subtracting terms needed to make the confining 
potential vanish at the origin: 



vr 



M 



C + M 



1 

+ - 



M 
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The remaining part of the effective Hamiltonian, V, contains, among other 
terms, emission and absorption of low energy gluons. Interactions which 
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change particle number enter at second order bound-state perturbation the- 
ory, which requires solutions to the nonperturbative three-body bound-state 
problem. 

V also contains a rotationally noninvariant part of the effective confining 
potential. Our choice of Hq does not lead to any first-order corrections to S 
states due to the rotationally noninvariant part of the potential, and for any 
I ^ state it minimizes the number of terms which give nonzero corrections, 
thus making the calculations easier. 

The spin-dependent part of the two-body effective operators is included 
in V, and it is treated in first-order bound-state perturbation theory. We will 
consider only the spin-spin hyperfine splitting in the ground state, because 
it can be calculated using the lowest order effective Hamiltonian [^J^. After 
a change of the spinor basis 0, the spin-spin part of the two-body effective 
interactions is: 



-'spin 



A{2MfaCFC^ 



6 TT 



4M2 



The function f{TZ) is rotationally noninvariant with respect to the angle 0. 
Its angular average is: 



ifm) 



sin7^ 2(1- cos 7e) 



This completes our discussion on Hq and V. 

We now want to solve nonperturbatively the eigenvalue problem for Hq: 

Hq\P)a = M^\P)a , 

where is the invariant mass of the bound state. We assume that the scale 
A is small enough so that the state is dominated by its qq component, i.e., 

d?K-^ dx 



\P)a 



i){K^,x)h^S\Q) . 



Let the mass of the bound state be 

= {2Mf + AME , 

which defines E. 

The eigenvalue problem for the Hamiltonian Hq leads to a dimensionless 
Schrodinger equation (^: 



d2 



1 



where 



^ + C -Fconf + T4oul(7^) 

dn^ \vr 



2maCF 



7/;(7e) = eV(7^) , 



2m{E - t) 
Z2 ■ 



In reference we show the dimensionless eigenvalue e for a few lowest 
lying states for c ranging up to 1. From the ratio of the splittings between IS 
and IP, and IP and 2S charmonium states, we find that c should be around 
0.6. This is the region where neither the Coulomb nor the confining potential 
dominate, and so there is not a simple Bohr scaling analysis which would 
assign powers of the coupling to momenta. For charmonium, we find values 

2 

of the quark mass Mc, the cutoff A and a = |^ so that: (i) c = 0.6 (or, 
equivalently, the ratio of the IS-IP splitting to the 1P-2S splitting is roughly 
correct), (ii) the mass of the ground state is Mis = 3.0 GeV, and, (iii) the 
mass of the lowest lying P state is Aiip = 3.5 GeV. Note that these values 
are reasonable approximations given the magnitude of known corrections. We 
obtain Mc = 1.5 GeV, a = 0.5 and A = 1.7 GeV. Similarly, for bottomonium 
we require c = 1.0, M.is = 9.4 GeV, and M.ip = 9.9 GeV leading to a bottom 
quark mass Mi, = 4.8 GeV, a = 0.4, and A = 3.5 GeV. It is important to note 
that this coupling need not run like a^jg or aiattice- 

With the parameters fixed, we can predict the hyperfine splitting in the 
charmonium ground state using eqn. (11). The function f{TZ) is rotationally 
noninvariant, but at c = 0.6 the violation of rotational symmetry is negligible. 
We predict that the splitting between the ground state vector and singlet in 
charmonium is 0.13 GeV, in reasonable agreement with experiment (0.118 it 
0.002 GeV). 

Next, we evaluate corrections due to the rotationally noninvariant part 
of the confining potential. There are first-order corrections to the P state, 
and the S state is corrected in second-order bound state perturbation theory. 
Corrections to the ground state are consistently a few percent even for c as 
small as 0.1. Corrections to the excited states at c = 0.6 and c = 1.0 are 
about 10%. This is a reasonable starting point because corrections of one 
higher power of a are of this same order. 

In conclusion, the logarithmic confining potential which arises at second 
order is a promising starting point for QCD calculations. Corrections to the 
energy levels due to rotational symmetry violating terms in this potential 
are negligible for the ground state, and for the lowest excited states they are 
small enough that corrections from higher order terms may restore rotational 
symmetry. This calculation is not intended to compete with phenomenological 
constituent quark model calculations. It is intended as an initial crude step 
toward and accurate first principles, light-front QCD calculation. 

Our work was supported by the National Science Foundation under grant 
PHY-9409042. 
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